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>D Abstract. Passive source localization is a well known inverse problem in which we convert 

" the observed measurements into information about the direction of arrivals. In this paper we 

^_^ focus on the optimal resolution of such problem. More precisely, we propose in this contribution 

Li^ to derive and analyze the Angular Resolution Limit (ARL) for the scenario of mixed Near- 

KH Field (NF) and Far-Field (FF) Sources. This scenario is relevant to some realistic situations. 

f/2 We base our analysis on the Smith's equation which involves the Cramer- Rao Bound (CRB). 

O This equation provides the theoretical ARL which is independent of a specific estimator. Our 

methodology is the following: first, we derive a closed-form expression of the CRB for the 

^vj considered problem. Using these expressions, we can rewrite the Smith's equation as a 4-th order 

►^ polynomial by assuming a small separation of the sources. Finally, we derive in closed-form the 

r~^ analytic ARL under or not the assumption of low noise variance. The obtained expression is 

l/"^ compact and can provide useful qualitative informations on the behavior of the ARL. 
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o 
^' 
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en 1. Introduction 

^^ Very few works are related to the study of the reahstic situation where there exists coexisting far- 

J> field (FF) and near-field (NF) sources [2] such as speaker localization using microphone arrays 

k>( and guidance (homing) systems. At the contrary, we can find a plethora of contributions on 

\^ the localization of far-field sources [1]. More recently, the problem of localization of near-field 

sources has been tackled in reference [H|5] for instance. In the context of the problem of source 
localization, one can see three contributions: (1) propose new efficient algorithms/estimators 
[1], (2) study the estimation performance independently of a specific algorithm thanks to the 
lower bound on the Mean Square error (MSE) [3l [5] and (3) derive and study the theoretical 
resolution, i.e., the minimal angular distance to resolve/discriminate two closely spaced emitted 
signals in terms of their direction of arrivals. Our contribution belongs to the third point. More 
precisely, based on the Smith's equation |6l [H |9] which involves the Cramer- Rao Bound (CRB) 
[3], we derive and analyze the Angular Resolution Limit (ARL) for the realistic scenario where 
we have two sources, one located in the far- field of the array and another in the near- field of the 
array. 



withs(t) = [si(t) S2{t)] 



T 



2. Model setup 

We consider some practical applications where the signals si(t) and S2{t) with t £ [1 : T] in 
which T is the number of snapshots received by an uniform linear array composed by L sensors 
are the mixture of near-field and far-field sources. More precisely, the (TL) x 1 observation 
vector is defined as y = y -|- e where e is the complex centered circular additive white Gaussian 

rAs(i) 

noise of variance o"^ and the (TL) x 1 noise-free signal is y = : 

[As(r) 

and A = [a(a;i) b(a;2i0)] where the signal sources Sm(i) are viewed as deterministic known 
signals and the steering vectors are defined by [a(a;i)]^ = e*"^^^ and [b(u;2, 4')]e = e'^i^^i+cpi ) -^vith 
i £ [0 : L — 1]. In addition, we assume wi ^ uj2. We define the determinist non-zero separation 
by (5 = UJ2 — uji- We focus our attention on the electric parameters but these parameters 
can be linked to the physical parameters, namely the DOA Om and the range r, according to 
Wm = —2TTd/Xsin{6m) and cp = 7rd?/{Xr) 008(62)'^ where d is the distance inter-sensor and A is 
the wavelength. 

3. Analytic ARL in the Smith's sense 

3.1. Analytical expression of the Cramer-Rao Bound 

The Cramer-Rao Bound (CRB) verifies the covariance inequality principle [3]. This bound is 
largely used in the signal processing community since it gives the best performance in term of 
Mean Square Error (MSE) at high Signal to Noise Ratio (SNR). Let cOm be an unbiased estimate 
of LOm, then 

E{{^m — <^m) } > [CJmm = [J ]mm (1) 

where E{.} is the mathematical expectation and J is the Fisher Information Matrix (FIM) 
defined by the Slepian-Bangs formula for a complex circular Gaussian observation y ~ 

CM {y, any. 

■^>« ^ >{(£)"^}^^|:«f<"''(S)''l^^<"} <^' 

where 1^ = [a(u;i) O] , || = [O b(w2»], |f = [O b(u;2,0)] where a^) (resp. b(^2,<A)) 
is the first-order derivative of a(a;i) (resp. h{u}2,(j))) w.r.t. the parameters wi, U2, respectively 
and h(uj2, (j)) is the first-order derivative w.r.t. the parameter (j). After some algebra, one obtains 
the CRB for the parameters of interest and the coupling term: 

CRB(^i) = ^|, (3) 

2L2||si||2-^^^ 

CRB(.;2) = y ^"^"^^^^ , (4) 

CRB(u;i,a;2) = -y ^ Q ' (^) 

where /? = ||s2||2 (^2 - g), Q = /3 (l2||si||2 - J^) - (c{5) - r?(5)g)', Lr = Y.tl P 



3.2. Derivation of the analytic ARL 

3.2.1. The Smith's equation: The Smith's methodology [B] provides the ARL as the solution 
of the following equation: 

CRB((5) = 52. (6) 

Expanding CRB((5) w.r.t. the localization parameters and the coupling term, we have [7]: 

CRB((5) = CRB(a;i) + CRB(w2) - 2CRB(wi, cja) (7) 

where the CRBs are given by expressions ([3|-([5]). 

3.2.2. Polynomial resolution of the linearized problem: Solving analytically equation ([6]) with 
M and ©-([s]) seems an intractable problem. Recalling that we assume that the separation 5 is 
small then a first-order Taylor expansions of C,{5) and r]{5) lead to the following approximations: 



r?(5) « ^{h{v{4>)+j5r{4>))} , C{S) « R{hiu{^) + j6v{m 



(8) 



in which n((/)) = Y2 



£=0 *- ^ 



-L-1 «3 



Eto^V-^' andr(0) = E 



L-l pA 



' . The above relation 



can be rewritten in a linear form w.r.t. the separation according to 

^{h{v{4>)+j5r{(t)))} = P-5Q, ^{h{u{4>)+j5v{(t)))} = P' -5Q' (9) 

^{hu{(j))} and Q' = —'^{hv{(j))}. Using these 



where P = 3f?{/ii;(0)}, Q = -9{/ir((^)}, P' = 
approximations, one can rewrite the CRB as 

2 



X 



(10) 



where, the Schur complement can be approximated as follow Q 



in which we 



>2(<5) pM 
[Pi{5) /3 
have introduced the two following polynomials P2{5) = aid + ao and P2{5) = 02^^ + ai5 + ao 



where 02 



_^ 



2112 ) '^i 



2PQ ^ 

6,2112 1 "0 



-C'2||si| 



L4||s2||2' "--L ~ L4||s2 

The linearized CRB expressions are now given by 



CRB(wi) 
CRB (1^2) ? 

CRB(a;i,W2)~ 



4^, ai = Q'-^^Q and ao = P'-gP- 



/3 



2 Q(,5)' 

(T^ i^2('5) 


2 Q(<5) ' 
ct2Pi(<5) 



(11) 



:i2i 



(13) 



2 Q(<5)' 
where Q[6) = P2{a) (3 - P^ (6) = {/3a2-al)S'^ + {/3ai-2aoai)5+l3ao-al. Consequently from (|6| 



and (11)-(13), the Smith's equation becomes 5'^ 



^ where Q'(5) = ^{f] + P^{5) + 2Pi{6)) 



del 



C20 + ci6 + Co with C2 = \a2, ci = \{ai + 2ai) and co = \{fi + ao + 2ao)- So, it is easy 
to see that the Smith's equation provides the ARL as the solution of a 4-th order polynomial 
according to 



R{x) = Q{x)x^ - 

where x is a free variable with qq 
where we assume j3a2 7^ a\. 



Q'{x) = x^ + g3X^ + g2x'^ + gix + go 



_C0_ 



/3a2-oJ 



T, gi 



Cl 



/3a2-aj 



7, 52 



l3ao-al-C2 
I3a2-a{ ' 



93 



(14) 

l3ai—2aoai 
/3a2-«i 



3.2.3. Analytic solutions of R{x): The resolution of a 4-th order polynomial is not 
straightforward but we propose a solution to this problem. More precisely, as CRB(5) = 
CRB(— (5), it is easy to see that if 5 is a solution of the Smith's equation then —6 is also a 
solution. This implies that if 5 is a root of R{x) then —5 is also a root. So, R[x) has four roots, 
namely {J, —5, ri, r2} and a decomposition of this polynomial into a product of monomial terms 
is given by 



R{x) = [x — 5){x + 5){x — ri){x — r2) 
= x^ — (ri + r2)x^ + (rir2 — 5'^)x'^ 
+ 5^{ri + r2)x -rir25^. 



(15) 
(16) 
(17) 



We can identify the coefficients of the polynomials (|14|) and (|17|) according to 



-(n +^2) 

rir2 — 5"^ 

(5^(ri +r2) 

— rir25^ 



I3a2-af ■ 
/3ao-a§-C2 
I3a2-a'j ■ 

/3a2-af' 
CQ 



(18) 



/3a2-af ■ 

Combining the second and the last equations, we have to solve the following polynomial: 



R'{x) = (/3a2 — ai)x + (/3ao 



"o 



C2)X 



cO. 



Polynomial R'{x) can be reformulated as a 2-rd order polynomial according to 



R'{z) = {I3a2 - a\)z^ + (/3ao 



"o 



The discriminant is 
A 



(^ao 



an 



C2 



1 + 4 



(/3a2 



02)2: — cO. 



a?) Co 



(/3ao 



an 



C2 



(19) 
(20) 

(21) 



o-2->0 



The study of the sign of the discriminant is not straightforward but observe that cq, C2 "— >" 0, 
so if the noise variance is low then the discriminant is given by 



A "'-^^ {(3ao 



"0 



C2Y > 0. 



(22) 



So, we know that if the noise variance is not too high, there exists two candidates for the 
ARL (or a double solution if the discriminant is zero) which are given by 



-(/?ao 



at\ 



C2)±VA 



2{/3a2 - ai) 



(23) 



where v A is given by the square root of expression (21). 



3.2.4. Analytic expression of the ARL: To discriminate the two possible solutions, we advocate 
that 6 



>o 



is a reasonable property. Thus, note thanks to (22 ), the sign in (23 ) must be chosen 

as positive to ensure 6 — J- 0. If not (in case of negative sign in (23)) the chosen solution will 
be meaningless. So, finally the ARL is given by 



-(/3ao -al -C2) + VA 
2(/3a2 - af) 



(24) 



To provide further simplifications, we can see that 

2(/3a2 - al)co 
/3ao - al-C2 



A « /3ao - Oq- C2 + 



(25) 



thanks to a first-order Taylor expansion of the square root \/l + x ~ 1 + l/2x for small x. Using 
the above approximation of vA, the ARL takes the simple following expression: 



Co 



/3ao - Oq - C2 



(26) 



So, we can see that 5 ^ 0{a). 



4. Numerical illustrations 

In this simulation part, we have considered an array constituted by L = 10 sensors. The sources 
are chosen to be close where Opp = tt/3 and 9nf = vr/3.1 and the range of this source belongs 
to the interval [0.62{d^{L - lf/T)X,2(f{L - if/X] |5j where d = 0.0125 m and the carrier 
frequency is /o = 10 Mhz (the wavelength is A = c/fo where c is celerity of the light). The 
modulus of the NF sources is fixed to ten times higher than the modulus of FF source. It is 
normal to assume that the source which is the closest has a higher power than the one which 
belongs to the FF. The number of snapshot is T = 100. On Fig. [T| we have reported the positive 
roots of each polynomials R{x), R'{x) and the analytic ARL given in expression (24). As we can 
see one root for polynomials R{x) and R'(x) are independent from the noise variance. These 
two roots have to be ignored. In addition, we can see that the one root for R{x) and R'{x) are 
equal and follow a decreasing function w.r.t. the inverse of the noise variance. Moreover, the 



analytic ARL given in expression ( 24 ) is in good agreement with these roots and thus assess the 



validity of the derivations given in the previous section. The ARL given in expression (24) and 



the approximated ARL under the assumption of low noise variance given in expression ( 26 ) are 



reported on Fig. 2. It is important to highlight the good accuracy of the proposed approximated 
ARL. 
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Figure 1. (a) Positive roots for polynomials R{x), R'{x) and analytic ARL given in expression 



(24) vs. the inverse of the noise variance, (b) ARL given in expression (24) and the approximated 
ARL under the assumption of low noise variance given in expression (26) vs. the inverse of the 
noise variance. 



5. Conclusion 

In this paper, we have derived and analyzed the Angular Resolution Limit (ARL) based on the 
resolution of the linearized Smith's equation for the new and realistic scenario where far-field 
and near-field sources are mixed. Generally speaking, the ARL is important and fundamental 
since this quantity gives the limit in the resolvability/separation of two closely spaced signals 
in term of their direction of arrivals. We show that for the chosen application, the resolution of 
the Smith's equation turns to be the selection of the "right" root of a 4-th order polynomial. 
This allows us to give a closed-form (analytic) expression of the ARL. 
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